arXiv:hep-ph/9808243vl 6 Aug 1998 


Applications of the Large Mass Expansion 


J. Fleischer Q A.V. Kotikov B0H and O. L. Veretin 00 


Fakultdt fiir Physik, Universitdt Bielefeld, D-33615 Bielefeld, Germany. 

The method of the large mass expansion (LME) is investigated for selfenergy 
and vertex functions in two-loop order. It has the technical advantage that in many 
cases the expansion coefficients can be expressed analytically. As long as only one 
non-zero external momentum squared, q^, is involved also the Taylor expansion 
(TE) w.r.t. small yields high precision results in a domain sufficient for most 
applications. In the case of only one non-zero mass M and only one external 
momentum squared, the expansion w.r.t. q^/M'^ is identical for the TE and the 
LME. In this case the combined techniques yield analytic expressions for many 
diagrams, which are quite easy to handle numerically. 


^ E-mail: fleischer@physik.uni-bielefeld.de 

^ Particle Physics Laboratory, Joint Institute for Nuclear Research, 141980, Dubna 
(Moscow Region), Russia 
® E-mail:kotikov@thsun 1.jinr.dubna.su 

Supported by Volkswagen-Stiftung under 1/71 293 
® E-mail:veretin@physik.uni-bielefeld.de 
® Supported by BMBF under 05 7BI92P(9) 



1 


1. Applications of the presented methods in particle physics 

Diagrams with one non-zero mass and only one external have wide 
applications in QED and QCD for both selfenergies and vertices. For ap¬ 
plications in the electroweak part of the fnll Standard Model (SM), due to 
the rather large spectrum of masses, in general one has to apply numer¬ 
ical methods. Often, however, approximations like small fermion masses 
mj/m| ~ 0, close vector boson masses (m^ — ~ 0 and large top 

mass m^/m^ ~ 0 etc. are considered. In these approximations the problem 
is often reduced to diagrams with one non-zero mass, which is a further 
reason for the interest in this type of diagrams. For applications in the 
electroweak part of the full SM, two calculations in the two-loop order are 
of particular interest: 

a) Z decay into bb (bottom, anti-bottom) 

The kinematics in this case is and Pi = P 2 = 0 (in the mi, = 0 

approximation). 



Fig. 1. Kinematics for the decay Z —> bb. 

We mention that the decay Z ^ bh was considered by a number of 
authors in the limit of large top mass, i.e. in the leading order only 
the top mass was kept non-zero and all other masses put to zero and even 
gi2 = = 0. This is an example of the above mentioned approxima¬ 

tions and here also the number of diagrams is quite low (of the order of 
10). Recently higher order terms in the large top mass expansion have 
been considered |§] (see also these proceedings). In Q an investigation of 
the precision of the higher orders in the large mass expansion has been 
performed for scalar diagrams occurring in this process and a repetition of 
higher orders in the large top mass for the full process is under consideration. 
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b) Anomalous magnetic moment (AMM) of the muon {g — 2) 

Here the kinematics is quite different, namely = 0 (actually q = 0) 
and Pi=P 2 = rnfi- In fact in this case the problem reduces to the calcula¬ 
tion of two-point functions. In both cases a) and b) we have one external 
variable only: and Pi = P 2 = rnli, respectively. 



Fig. 2. Kinematics for anomalous moment {g — 2)^. 

In the case of the AMM our first step was dedicated to an automatic 
set up of the contributing Feynman diagrams. For this purpose the package 
TLAMM was developed Q. Also the large mass expansion (LME) ||5| ap¬ 
plied to a “toy model” for the {g — 2) was described in |Q], i.e. for diagrams 
of the selfenergy type. In this model 40 diagrams were contributing to the 
two-loop AMM of the muon. 

In general, however, in the above cases many more, i.e. of the order 
of 1000 diagrams contribute and naturally it is necessary to investigate 
diagrams with various non-zero masses. 

As was already started with TLAMM, in the full SM it is even more 
required to produce the Feynman diagrams automatically. Instead of 40 
diagrams, e.g., in the above mentioned “toy model” in the full {g — 2)^ cal¬ 
culation 1832 diagrams contribute. Such a more general program, called 
DIANA (for Diagram ANAlyser) Q, written in C and making use of 
Nogueira’s QGRAF Q, was also developed in our group. It produces FORM 

input for the Feynman diagrams according to the Feynman rules for fur¬ 
ther evaluation. 
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2. Evaluation of scalar diagrams 

In the following we shortly review the main techniques we applied, 
namely the Taylor expansion (TE) and the LME [^. The efficiency 
of both approaches will be compared. 


2.1. The Taylor expansion method 


Taylor series expansions in terms of one external momentum squared, 
say, were considered in 0, Fade approximants were introduced in |jl^ 
and in Ref. Q it was demonstrated that this approach can be used to 
calculate Feynman diagrams on their cut by analytic continuation. In the 
case of a three-point function like Z —> 56 we have two independent external 
momenta in d = 4 — 2e dimensions. The general expansion of (any loop) 
scalar 3-point function with its momentum space representation C{pi,p 2 ) 
can be written as 


C{P1,P2) = £ aimn{piy{pl)"^{PlP2y 


( 1 ) 


l,m,n=0 


where the coefficients aimn are to be determined from the given diagram. 
Considering Z bb with mb = 0, i.e. Pi = P 2 = 0; which is a good 
approximation, only the coefficients ooon are needed. For the calculation 
of the Taylor coefficients in general various procedures have been proposed 
|Q,|12],||^. These methods are well suited for programming in terms of a 
formulae manipulating language like FORM. Such programs, however, yield 
acceptable analytic results only in cases when not too many parameters (like 
masses) enter the problem. Otherwise numerical methods are needed |]l3, 

0 . 


In the case of only one non-zero mass and only one external momentum 
squared, indeed the case with the least nontrivial parameters, for many 
diagrams analytic expressions for the Taylor coefficients can be obtained. 
In Sect. 5 we present some recent results. 

For the purpose of calculating Feynman diagrams in the kinematical 
domain of interest it is necessary to calculate them from the Taylor series 
on their cut. This is performed by analytic continuation in terms of a 
mapping. 

Assume, the following Taylor expansion of a scalar diagram or a particu¬ 
lar amplitude is given C'(pi,p2) • • •) = Z)m=o = f{y) and the function 

on the r.h.s. has a cut for y > yo- 

The method of evaluation of the original series consists in a hrst step in 
a conformal mapping of the cut plane into the unit circle and secondly the 
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reexpansion of the function under consideration into a power series w.r.t. 
the new conformal variable. A variable often used is 


1- yi- y/yp 
1 + \/i - y/yo 






Fig. 3. Conformal mapping of the complex y-plane into the w-plane. 


By this conformal transformation, the y-plane, cut from i/q to +oo, is 
mapped into the unit circle (see FigJU) and the cut itself is mapped on its 
boundary, the upper semicircle corresponding to the upper side of the cut. 
The origin goes into the point cu = 0. 


After conformal transformation it is suggestive to improve the conver¬ 
gence of the new series w.r.t. a; by applying the Fade method |16|,[17|. 
A convenient technique for the evaluation of Fade approximations is the 
e-algorithm of [1T| which allows one to evaluate the Fade approximants 
recursively. 


2.2. Large mass expansion (LME) 

In particular for the evaluation of diagrams with several different masses, 
one of which being large (like the top mass, e.g., we use the general 
method of asymptotic expansion in large masses i- For a given scalar 
graph G the expansion in large mass is given by the formula 

Fciq, M, m, e) ^ Fc/^iq, m, e) o Tqj^^-yF.y{q^, M, mG, e), (3) 

7 

where y’s are subgraphs involved in the LME, G/y denotes shrinking of 
7 to a point; Fly is the Feynman integral corresponding to 7; Tq^^rn^ is 
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the Taylor operator expanding the integrand in small masses {m-y} and 
external momenta {( 77 } of the snbgraph 7 ; o stands for the convolution 
of the subgraph expansion with the integrand Fqj^. The sum goes over 
all subgraphs 7 which (a) contain all lines with large masses, and (b) are 
one-particle irreducible w.r.t. light lines. 

For the Z ^ bb decay we have for the on-shell Z’s. Fig.4 shows 

diagrams with two different masses on virtual lines, one of which is a top. 
W and Z are the gauge bosons with masses M\y and Mz, respectively; cj) is 
the charged would-be Goldstone boson (we use the Feynman gauge); t and 
b are the t- and b-quarks. 



(p2-pl) 


(p2-pl) 


(p2-pl) 


case 7 case 7.1 case 7.2 


Fig. 4. Two-loop diagrams with two different masses in internal lines arising in 
the process Z bb. The notation for the diagrams is chosen according to ||]. 




case 7 


case 7.1 

V V V 


case 7.2 



Fig. 5. The structure of the LME, see explanations in the text. 


Applying the method of the LME, Mp/ and Mz are considered as small. 
In the framework of this expansion contributions from additional subgraphs 
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are to be taken into account together with the Taylor expansion of the ini¬ 
tial diagrams w.r.t. external momenta and light masses. These are shown 
in Fig.5. Bold, thin and dashed lines correspond to heavy-mass, light-mass, 
and massless propagators, respectively. Dotted lines indicate the lines omit¬ 
ted in the original graph T to yield the subgraph 7 . T/y (see (^)) then 
consists out of all the dotted lines after schrinking 7 to a point. 

These subgraphs restore the analytic properties of the initial diagrams 
(like logarithmic behaviour). In other words, the Taylor expansion of the 
initial diagrams produces extra infrared singularities which are compensated 
by singularities of the additional subgraphs so that only the singularities of 
the original diagrams survive. 

At this point it also becomes clear what the difference is between the 
small-g^ expansion and what is called here the LME; in the former case 
we assume all masses large, i.e. (f -C while in the latter 

case only mt is considered as large and all other parameters small, i.e. 

<C m^. In this sense both methods are LM expansions. The 
technical advantage of the second method is, however, that only bubbles 
with one mass occur, which can be expressed in terms of T functions, while 
in the case of the small expansion bubbles with different masses are 
involved, which are much more difficult to evaluate, in fact only numerically. 
Of course also the number and structure of the subgraphs is different in the 
two cases. 

Finally the LME of the above diagrams has the following general form: 


= — 



Til ( 4 ) 


i-\-j = n 


where m is the highest degree of divergence (ultraviolet, infrared, collinear) 
in the various contributions to the LME (m < 3 in the cases considered). 

and /mf are considered as small parameters. are in general 

complicated functions of the arguments, i.e. they may contain logarithms 
and higher polylogarithms. 

In the next Section we present results and compare the two methods 
described here. 


3. Results 

The small momentum expansion of cases 7, 7.1 and 7.2 is described in 
detail in Ref. [^]. The additional subgraph arising in these cases is the 
same irrespectively of the mass distribution and is shown in Fig.5 {case 7). 
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Fig. 6. Results of the LME for the real finite parts of diagrams 7.1 and 7.2 . Solid 
curves represent the small-g^ expansions, middle-dashed the leading term of the 
LME, short-dashed the sum of 11 terms in the LME, long-dashed the [5/5] Fade 
approximant from the LME. 



q^lra\ 



In the case of the LME two additional subgraphs arise in each of the 
cases 7.1 and 7.2. Beyond that these sets of additional subgraphs are also 
of quite different structure. Furthermore the 2"^^ and 4**^ graphs in the 
row (see Fig.5) produce 1/e^ terms which cancel, however. Since there are 
no UV divergences, these must be mixtures of infrared and collinear ones. 
They determine the highest degree of divergence in these cases and thus the 
highest power of the logarithm as discussed in (4). 

Our numerical results for cases 7.1 and 7.2 are presented in Fig.6, and 
for in Table 1. In the figures we show the small g^-expansion 

(solid line) in comparison with the lowest order approximation (middle- 
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dashed) and the sum of terms (small-dashed) with N=9, see (4). The scale 
parameter /i = mt, i.e. only k=0 contributes in (4). We see that up to 
the sum of 11 terms agrees quite well with the small g'^-expansion, 
while for higher the agreement very quickly worsens. For this reason 
we formally apply also the Fade summation technique. With 11 terms in 
the series, a [5/5]-approximant (long-dashed) can be constructed. It is seen 
that indeed this improves the situation considerably up to g^ = 4M| though 
for a better agreement many more terms in the LME would be needed. In 
the small g^-expansion only 9 terms were taken into account, i.e. a [4/4] 
approximant is calculated. From Table 1 we see that for g^ = indeed 
a rather precise result can be obtained with 11 terms from the LME, in 
particular if Fade’s are applied [^. 

^(0) Table 1 corresponds to the lowest order of the diagrams (N=-l 
or 0, respectively) and to the order of expansion performed in Ref. |2|. 


Table 1. Values of diagrams at = 90^, /r = toj = 180, mw = 80. 



case 7.1 

case 7.2 

p{^) 

6.4 

-i30.0 

9.2 

-i24.5 

pW 

11.1 

-il7.93 

9.4 

-*44.848 

piw) 

10.1 

-il7.95218 

10.122 

-*44.84523 

[5/5] 

9.996 

-il7.9528 

10.189 

-*44.842 

small-g 

9.992668259 

-il7.95215366 

10.193902102 

-*44.845273975 


4. Bremsstrahlung 

Of particular interest is the calculation of Bremsstrahlung in terms of 
the LME. Let us consider again the decay Z —> bb. The kinematics of the 
(gluon-) Bremsstrahlung for this process is given by 

q^Pi-P2+P3, pI=pI=pI = 0- 
Thus we have 3 invariants 

PlP2, P2P3, P3Pl- 

We are interested in the integrated bremsstrahlung. The phase space 
measure in d dimensions can be written as 

1 2(27r)''-2g2 /g2y-" 

(27r)3<i-3 Y{d - 2) 32 2 ) 

[ dxdyx‘^~^{l - - yfl'^~'^ 

Jo 


( 5 ) 
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and the invariants can be expressed in terms of x, y as 

PiP 2 = -jx{l - y ), 

P2P3 = yil - x ), 

pwi = -j^y- 

For the 1-loop amplitude we apply the LME w.r.t the top mass . For 
the diagrams with Mw and mt mass in virtual lines the desired expansion 
looks like 



where fijkn ai'e cubic polynoms of \og{rnl/y?) with coefficients being func¬ 
tions of i.e. 


fijkn = a log'^ + blog^ -^ + clog^ + d. 

y y y 

After expansion the integration over the phase space is trivial (it can be 
done completely in terms of F-functions). 


5. Diagrams with only one non-zero mass 

So far the method of expansion is considered as semianalytic in the sence 
that only a limited number of coefficients can be obtained explicitely. In 
this Section, however, we want to go one step further. We calculate the 
first few coefficients of the expansion of a diagram with only one non-zero 
mass by means of the LME |^] method. The method of differential equa¬ 
tions 1 ^ then yields an idea, like in 0 , what the general analytic form 
of these coefficients might be, providing some “basis” in terms of which 
they might be expressed. The Ansatz of equating the explicit coefficients 
obtained from the LME to a linear combination in terms of the basis ele¬ 
ments, yields a system of linear equations, which can be solved to yield the 
desired representation of the coefficients. 

The main problem in this approach is the choice of the basis elements. 
We start from so-called harmonic sums which are particularly relevant for 
moments of structure functions in QCD (e.g. |^^). These functions are 
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directly related to (generalized) poly logarithms |21] and therefore it is not 
surprising that they appear in the analysis of massive diagrams. However, 
not every massive diagram even of self-energy type can be expressed in terms 
of these sums. Apart from harmonic sums we introduce a new type of sums 
which we call W- and H-sums. 

Our starting point is the LME of the diagrams. For the particular 
diagrams under consideration it was discussed in detail in Q. Here we just 
note that the result of a large mass expansion for these diagrams reads 
(z = jrv?) 


•^ = i E E 1 1: . 

' n>l j=0^ 


(7) 


fc =0 


where a is the dimensionality of the diagram, iz and v independent of n 
are the highest degree of divergence and the highest power of log{—q‘^/m‘^), 
respectively (in our cases a;,z^ < 4). The coefficients j ^ are of the form 
ri + C2'>'2 + ... + with Tc being rational numbers and Cc = C(c) is the 
Riemann ('-function. 

Series (j^ always has a nonzero radius of convergence, which is defined 
by the position of the nearest nonzero threshold in the g^-channel. For 
brevity we shall call m-cuts (2m-cuts, etc.) possible cuts of a diagram in 
corresponding to 1 (2, etc.) massive particles in the intermediate state. 

We start from diagrams with the simplest threshold structure, i.e. with 
m-cuts and 0-cuts. It turns out that all of these with only one massive 
line are expressible in terms of harmonic sums Sk{n — 1) = 1/j^ a-nd 

alternating harmonic sums Ka{n — 1) = (“)'^'*'^/j'^- 

As a first example consider the 2-loop 2-point function of Fig. 7, J(a), 
which was considered in |22]. Using the standard large mass expansion 


technique one can get the first few coefficients of the expansion of this 
diagram in powers of z = jrr? 


1 


^ n=l 


^ OnZ"- = 


-K [ 2 ( 2 ^; + ((2 + - 


2 :^ -b 


2 ^ 1 


Z'* -b . . . 


( 8 ) 


where (2 = C(2) is the Riemann (-function. 

Our first step is to find an expression for the higher order terms of the 
series (P) as a functions of n i.e. = a{n). 

To achieve this let us search for as a linear combination of harmonic 
sums. The rule is that one has to take into account all possible products of 
the type CaSbi ■ ■ ■ Sb^./n^ with the ’transcendentality level’ (TL) a-b6i-b.. .-b 
bk + c = 3. It is obvious that one can exclude (3 from the very beginning 
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Fig. 7. Three diagrams, selfenergy and vertices, for the TE of which results are 
reported. 


since it never appears on the r.h.s. of (g). Thus we have the following 
Ansatz for 


an = —Xx+C,2SlX2+S^X2,+S2SiX4^+SlX^+‘- 
n 


where xi,...,xg are rational numbers independent on n. The argument 
of the S'-functions is n — 1. We refer to the structures in as ‘basis ele¬ 
ments’. Indeed, the functions CaSbiin — 1)... Sb^.{n — l)/n‘^ are algebraically 
independent. 

Inserting the expression for Un from into the l.h.s. of @ and equating 
equal powers of z, we obtain a system of linear equations for the Xi. An ex¬ 
plicit computation ensures that the system can be solved in terms of rational 
numbers for the Xi The solution is xi = 2,x% = 2,X8 = —2, X 2 , 3 , 4 , 5 , 7,8 = 0 
i.e. the answer for the diagram at hand is 


J{a) 


oo 


E- 

n=l 



n 


2h(!^y 

I 


( 10 ) 


Equation (P) can be considered as expansion of the general coefficient 
an in terms of a ‘basis” with (unknown, rational) coefficients Xj. Now the 
question arises: what are the basis elements in general? So far we are lacking 
rules for predicting a basis of a given diagram. The power of the method, 
however, is that given a set of basis elements for one diagram, it can be 
used to find the solution for other diagrams. Often, though, one has to 
‘generalize’ already known (harmonic) sums in a basis. 

As an illustration consider the diagram shown in Fig. 7, Its large 
mass expansion looks like 


J{b) - 


(q 


^ f; f r(2) log2 (-^) + rW log{-z) + , 


n=l 


( 11 ) 
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with the r’s being rational numbers. It is obvious that one can search for 
a solution for each of the r’s independently. Again we use the same set of 
functions (l/n“ and Sb) as above but now with different transcendentality 
level(s). The system of equations has a solution only if we add the factor 
(—)”■ which can be seen easily by inspection of the series. At the end we 
arrive at 




n=l 


-Z 




n 


n 


n 


n 






n'^ 


( 12 ) 


where as above we take all harmonic sums [Si) with the argument n — 1, 
which is omitted. 

We point out that in the part without log(—z) has basis elements 
QaSb/n^ obeying a + 6 + c = 4 (i.e. it has a basis with TL= 4). Terms 
proportional to log(—z) are of 3rd level while those proportional to log^(—z) 
of the 2nd. One can say that log(—z) itself is of 1st level and each log“(—z) 
reduces the level of basis elements by a units. This is the general behaviour 
for all diagrams we consider. The same rule applies to UV and IR poles 
1 /e if they are present in a diagram i.e. each l/e“ reduces the level of basis 
elements by a units. 

We also stress the presence of the factor (—)”■. Strictly speaking the 
basis now is not CaSb/n^ but rather {—)'^CaSb/n^- 

Finally we observe the following recipe of generalization,i.e. elements of 
higher TL (a + h) can be constructed in the form: 


n , 

Sa,b{n) = -;^Sb{j - 1), (13) 

j=i J 

(-V 

Ka,b{n) = (14) 

7 = 1 J 

(15) 

The structure of the lowest level for diagrams having both m- and 2m- 
cuts is more complicated. We find it from the nonplanar graph shown in 
Fig. 7, J(c). Namely this diagram has a double collinear pole and we can 
easily perform the integrations in calculating the 1/e^ contribution with the 
result 


1/£2 

J(c) = - 


dcxi do.2 


£2(q,2)2 _ a2)z)(l - (1 - ai)a2z) 
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E 


(n !)2 


( 2 ^- 1 )' UU'-yj 


e2(g212 


E 


^ /2n\ ^g H^i(n-l) 


n=l 


n 


n 


(16) 


If we assign to the factor (^”) ^ the 0th level, then the expression on the 
r.h.s of ([T^ ) has the correct 2nd level as it should be for the double pole 
part of a vertex function. The generalization is done according to 


= Er^)4*560'-1)- 




The result for the diagram under consideration finally is 


(17) 

(18) 


j(c) — 


(q 


2\2 


E^ 

n=l 


2n\ 1 


n n 


-^lVi + 


1 r 


e L 


- 4 IT 2 - 121Ti,i - 125ilTi - 1652 
— 18(^2 bid — 4IT3 — 853 T 121bi,2 — — I2IT3 

-165i,2 - 85i 1T2 - 245ilTi,i - 165i52 - 125?lTi I, (19) 


where 


Sa,{b+c){n) = - ^)Sc{j - 1 ) , 


j=l 3 


Wa,(h+c){n) = ^ ~ 


( 20 ) 

( 21 ) 


all have TL= a + b + c. 

For completeness we mention a further type of basis elements which we 
discovered from the application of the differential equation method fTT 

"-"■sBJf 


( 22 ) 
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Further basis elements were obtained in a similar fashion. The result for 
a great variety of diagrams will be published elsewhere. 


6. Conclusion 

The purpose of many of our investigations is to test various methods for 
the calculation of two-loop diagrams needed for the evaluation of measured 
processes in the full SM of the electroweak theory. Comparing different 
methods yields information about their precision and possible applicability. 
In particular we have forced the developement of packages for the large 
mass expansion and found that in the domain of interest for the processes 
under consideration this method can be used with reasonable liability. For 
the case of only one non-zero mass a completely new approach was found, 
namely “guess and verify”, to set up analytic expressions for many diagrams 
of interest. 
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